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Abstract: We calculate the twist-4 corrections to the integral of gi{x, in 
the framework of QGD sum rules using an interpolating nucleon held which 
contains explicitly a gluonic degree of freedom. This information can be 
used together with previous calculations of the twist-3 contribution to the 
second moment of g 2 {x) to estimate the higher-twist corrections to the Ellis- 
Jaffe and Bjorken sum rules. We get f^‘^\proton) = —0.037 ± 0.006 and 
(neutron) = —0.013 ± 0.006. Numerically our results roughly agree with 
those obtained by Balitsky, Braun and Kolesnichenko based on a sum rule 
for a simpler current. Our calculations are far more stable as tested within 
the sum rule approach but are more sensitive to less well known condensates. 
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Spin effects in strong interaction high energy processes are one of the best 
tools to study QCD both in the perturbative and non perturbative regime. It 
has become clear that the reliable determination of the dependence, both 
due to radiative corrections and due to higher twist contributions is a central 
task of QCD theory. The dependence of spin variables is in general more 
benign than for unpolarized quantities (e. g. the anomalous dimension for 
the Bjorken sum rule vanishes), allowing to extract very interesting infor¬ 
mation from data taken at as low as 1 GeV^. While the EMC [jl[ and 
SMC [|] experiments have still comparatively large , SLAC [§ data are 
taken down at a rather low mean of about 2 GeV^. The proton data 
indicated a disagreement with the polarized-proton sum rule, the Ellis-Jaffe 
sum rule |^, leading to a lot of excitement in the high-energy physics com¬ 
munity. Since experiments now not only provide data for the spin-dependent 
structure function of the proton gi{x) but also of the neutron gi{x) one of 
the most solid predictions of QCD, the Bjorken sum rule 0, is tested experi¬ 
mentally. However, for low this sum rule (strictly valid in the asymptotic 
Bjorken limit) receives corrections. The most familiar are those from per¬ 
turbative QCD, calculated for the leading twist term up to order ag Q and 
estimated to order ag with an rough estimate of order a| 0. Higher-twist 
corrections are given in terms of complicated hadronic matrix-elements and 
are suppressed by powers of Q^. Since experiments still are not in a save 
region of asymptotically large all this corrections have to be examined 
carefully to give a complete picture of the dependence of the Bjorken sum 
rule 0. 

The first moment of gi{x, Q^) at hxed is given by 0 

£ dx Q^) = bl») + ^ pp) + 4d<^) + 4/P)) + O(^) , (1) 

In the above formula we have not included higher order corrections to the 
coefficient functions which can be written as power series in ag{Q‘^) [| In the 
present paper we focus on the correction proportional to which is dehned 
by the matrix element of the twist-4 operator 

_ 0^= g(0)^7G./3(0)7^(0) (2) 

^The radiative corrections to the leading-twist singlet part were calculated to order a| 
0 and estimated to order Og |ll|] . 
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in the nucleon state IPS') of momentum p and spin S', = — S'^ = 

p - ^ = 0, 

(pS'|O<,(0)|pS') = . (3) 

Here and in the following we temporarily neglect the normalization-point 
dependence of the operator 0„. 

The reduced matrix elements dP''^ and can be expressed through the 
second moments of the polarized nucleon structure functions pi(x) and g 2 {x) 

W%- 


/ 


dx x‘^g 2 {x) 


J dx x‘^gi{x) 



(4) 


While may be taken directly from experiments, the matrix elements 
and were hrst estimated by Balitsky, Braun and Kolesnichenko JT^ 


using the QCD sum rules techinque. In our previous work [T^ we presented 
an independent sum rule calculation of which essentially conhrmed the 
values obtained by BBK. The calculation of presented in this paper 
completes therefore the program of estimating the leading power corrections 
to the sum rules for / gpx, Q^)dx. 

Before we are going to dwell on the details of the calculation we would like 
to stress that the leading higher-twist matrix elements describe fundamental 
properties of the nucleon. The twist-4 operator eq. is a measure for the 
contribution of the collective gluonic held to the spin of the nucleon. Writing 
the dual held strength tensor in its components we get 


(pS| - B\j\ + (]a X = 2m5,/P>S" 


( 8 ) 


where the quark-current is denoted as = —gq'j^Pq and and are the 
colour magnetic and colour electric helds. In the rest system of the nucleon 
an analogous relation holds for the twist-3 operator which determines d*-^^ 

+ iu X EaTIpS) = . (6) 

Knowledge of and then allows to estimate magnetic and electric held 
contributions to the spin separately. 
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In the usual approach nucleon matrix elements of local operators can be 
extracted from a three-point correlation function 

nr(p) = J d^y{0\T {r]{x)Or{y)r]{0)} |0) (7) 

which involves an interpolating current ri{x) with a certain overlap A between 
the state created from the vacuum by ri{x) and the nucleon state 

(0|?7(a;))|pS') = Xu{p, S) exp (ipx). (8) 

The overlap integral can be determined from an additional two-point corre¬ 
lation function 

U{p) = i J d^xe^P^{0\T {p{x)r]{0)} |0). (9) 


In practical application it is often advantageous to consider the ratio of three- 
and two-point correlation functions such that the A-dependence cancels out. 

For QCD sum rule calculations of nucleon properties the standard choice 
for T]{x) has been for a long time the three-quark current introduced by Ioffe 



pi{x) 


u'^{x)C'y^u\x) '^^Yd‘'{x)e' 


abc 


( 10 ) 


which was used in the calculation of and by Balitsky, Braun and 
Kolesnichenko (BBK) [^. As explained in [^, for the investigation of 
operators which, like eq. (|^), contain explicitly gluonic degrees of freedom 
it is very useful to match these by an interpolating nucleon current that 
contains gluonic degrees of freedom as well. For such a proton current we 
chose 



( 11 ) 


where 

and 


[G“^{x)d(x) 

ri'Six) = £“‘' (u'‘(x)C',^d\x)) 757'‘<£«£ [G"'’(l)u(a:)‘ 


( 12 ) 

(13) 


This current was hrst studied in and tested in the calculation of the nu¬ 


cleon gluonic form factor and the total momentum fraction carried by gluons. 
Next the twist-3 correction occurring in the expansion of / gi{x)dx was 
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predicted |^. Note that the current (pN]) may be regarded as the leading ex¬ 
pansion term of a non-local version of the classical three-quark current ( 0 ). 
Working with non-local sources has become also popular in lattice-gauge the¬ 
ories where stability can be increased by taking a number of derivatives of 
the quark-helds. 

Using lattice-gauge theory the three- and two-point correlators and 
(H) can be calculated directly from hrst principles. Such a project is pursued 

a 


by the Jiihch-group 


The basic idea of the QCD sum rule technique on 


the other hand is to employ the duality between the hadronic and partonic 
representation of a correlation function and to extract the quantity of interest 
by demanding that both descriptions match each other at some intermediate 
scale. If —p^ is sufficiently large the main contribution comes from small 
a:-distances of order of ~ l/(—p^). The consideration of contributions 
from different t-channel distances is more involved. When y‘^ < l/(—p^) the 
standard machinery of the short-distance expansion is applicable resulting 
in the known expansion in terms of quark and gluon condensates. However, 
the contribution from large has to be accounted for separately. 

The solution to this problem was first formulated by Balitsky [j^. The 
Operator Product Expansion (OPE) of a three-point correlation function has 
a twofold structure. Terms of the first type come from the region ~ ~ 

l/(—p^) and are proportional to vacuum expectation values (VEV) of local 
gauge-invariant operators multiplied by coefficient functions depending on p^. 
In the following we shall refer to these terms as to local power corrections 
(LPC). 

Terms of the second type called bilocal power corrections (BPC) originate 
from distances ~ l/(—p^). To treat such contributions properly one 

should expand the time-ordered product of nucleon interpolating currents 


r(,„(a:)^G(0))=i:C“(l)O„(0), 


(14) 


in a series of local, gauge-invariant operators O„(0) of increasing dimension n. 
When this expansion is inserted back in (|^) it results, together with standard 
LPC, in the following general form of the OPE of the three-point correlator 
( 0 ) 

nr(p) = I] C^^ni.p) <Ol>+Y. Cn^(p)nF,n(0) • (15) 


5 










The bilocal power corrections are determined by the long-distance, non- 


perturbative contributions to correlation functions at zero momentum |19 


n 


BL 

r,n 


(«)=*/ d‘ye'^{0\T {Or(!/)O„(0)} |0>, 


(16) 


with g = 0. Special care has to be taken to properly eliminate the pertur¬ 
bative short-distance singularities from these correlators which are already 
included in LPC coefficients Cp„(p). We stress again that the essence of the 
expansion (^) is the separation of different scales: vacuum expectation val¬ 
ues of local operators {O^) and correlators describe long distance effects 
while the coefficients and receive contribution only from highly vir¬ 
tual quark and gluon fields which propagate for small distances. The OPE 
of the correlation function Ifp must be written as the sum of both LPC’s 
and BPC’s, as above. We note that in general only this sum has a physical 
meaning and is independent of the regularisation scheme. Contrary to the 
case of our previous calculation of the twist-3 correction fl^, the BPC’s do 
play a crucial role in the analysis of the twist-4 matrix element. 

The right hand side of (J^) can be decomposed into different Lorentz struc¬ 
tures according to 

i J J A(0|T {vi{x)q{y)gG^f,{y)'y^q{y)ric{0)} |0) 

= + • • • , (17) 


where the ellipses stand for other terms which can be eliminated taking the 
trace of ([l7|) with an appropriate projector. The invariant functions W*, 
i = A,B can be represented as spectral integrals 


= / ds -, 

J s — 

+ subtractions. 


aA'{s — m%) -|- — m%) + 0(s — so)p‘^(s) 


(18) 


We have accepted the conventional “resonance plus continuum” model of 
the spectral density with the continuum density p^(s) dual to all graphs with 
non-vanishing imaginary part at —p^ —>■ oo. The constant a* which stands in 
front of the double-pole term is proportional to while the single-pole term 
is determined by nucleon-to-continuum transitions and has to be eliminated 
from the final answer. In principle the information about the magnitude of 
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can be extracted either from W'^ or from W^. In practical calculations it 
is advantageous to consider IV-^ because of its lower dimensionality and this 
structure was chosen for the analysis presented in []^. We realize that the 
QCD sum rule approach involves a number of approximations the accuracy 
of which is sometimes difficult to assess a’priori. To get a better feeling of 
the intrinsic uncertainties of the whole method we have decided to analyze 
which leads to a sum rule which is more sensitive to the higher mass 
region in the spectral representation (|T^. The constant coefficient if front of 
the double-pole nucleon contribution to can be found to be equal to 

= -2f^^'^\iXGm% . (19) 

Here, A/ and Xq are the overlap of the Ioffe current (|I0|) respectively the 
quark - gluon current (0) with the nucleon, while rriN denotes the nucleon 
mass. 

The invariant function can be easily projected out of the three-point 
correlator (|^) by taking the trace with and choosing momentum 

such that p ■ S = 0. The expansion of in LPC is standard and we are 
not going to dwell on the details here. The net result can be written as the 
following formula: 

4VkL^Pc(P^) = 

dl^(-p^)^log(-pV/^^) + B— < qq 

71° TT 

+cX < !</ggg > (-f) iog{-PlP) 

+F—rn^Q < qq 

^ gg >^< g'^GG > ^ qq m* 

_p2 _p2 

where the numerical coefficients are given in Table |^. Note that the coefficient 
in front of the gluon condensate (GG) vanishes. 

The calculation of BPC’s is more involved and will be described next 
in some details. Let us consider operators which may contribute two-point 
correlators (jl^) of dim-6 and dim-10 to the sum rule. The expansion of 
T{rij{x)'i]G{0)) in local operators leads to the following series: 


(-p2) log(-pV/^^) 


( 20 ) 
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71 

B 

G 

D 

G 

H 

d 

-1/2592 

-4/3 

0 

88/81 

5/81 

1/9 

u 

-1/3240 

-20/27 

-23/2304 

1411/2592 

1/81 

-1/54 


Table 1: Numerical coefficients corresponding to LPC contributing to the sum 
rule eq. . The upper line gives the values for the twist-4 operator involving d 
quarks, the lower line for the corresponding operator with u-quarks 


+ ^^^75 [* ( 0 ) - ( 0 ) 

+ml <qq> ^^7¥>75 2 (0) - ^ (0) 

■' {d^^D^d) ( 0 )' 

' 27^ (5a;‘^a;a - 7v.75 (m7“75m) (0) 

' nn ^2 rn'l 


2 ^2x^Xo, 

-- ' ' 27'kx^ 

< qq > Jl^^*^5 [(ugG"^D“u) (0) + 2 {dgG^^D“d) (0)] + 


The ellipses on the right hand side denote gamma-structures other than 7(^75 
that do not contribute to . Furthermore in the above formula we did not 
write explicitly the operator < qq > {qa^’^gG^'^D°‘q){0) due to the compli¬ 
cated and lengthy structure of its coefficient. In the end its contribution to 
the hnal sum rule turns out to be small. 

The relevant BPC arise from the correlation function of the local oper¬ 
ators present on the right-hand side of eq. ( |^) with the twist-4 operator 
which dehnes at zero momentum. In the particular case of operators 
which are proportional to equations of motion of QCD, like the operators 
q'y^'J^Ggaq, q'y 5 (J^xDaq and q'y^Daq, the correlation functions can be evalu¬ 
ated by means of exact low-energy theorems. Using the functional integral 





























representation it is easy to obtain the following Ward identities, see e. g. 


i 

i 

i 


d%{0\TO^iy) (qj^Daq) ( 0 )| 0 ) = 0 

d^y{0\TOa{y) {q^i^cr^pxDaq) (0)|0) =ml<qq> ^ {gaxga<t> 

47r 

d%{0\TO„{y) {q'j 5 YgGgaq) ( 0 )| 0 ) = —as < qq gaa 


gcripgaX^ 


( 22 ) 


To avoid misunderstanding we note that the above identities should be un¬ 
derstood in the following way. The large-distance contribution to the Ihs is 
equal to the large distance contribution to the rhs which is just determined 
by a non-perturbative VEV of the corresponding operator. 

The other remaining two-point correlation functions at zero-momentum 
transfer cannot be evaluated exactly. Instead, they may be estimated by 
considering additional two-point sum rules. Since ultimately one is interested 
in their long-distance behaviour, we focus on the contribution coming from 
the lowest possible intermediate state - the massless, chiral, pseudoscalar 
meson. One can dehne the overlap of the operators given in (|^) with a pion 
state as: 


(0|g7/.75gk(p)) = iUp^ 

(0|g^G'^^7^g|7r(p)) = ifJlp,, 

{^\qgGi,yD''q\TT{p)) = ifJlPf, 

{0\q(^i,egG^‘"D„q\n{p)) = ifJlp^ ( 23 ) 


The constants d-,,, and d-^ can be calculated by evaluating two-point cor¬ 
relation functions in the sum rule framework at non-zero momentum. We 
stress that to avoid double counting special care has to be taken to sub¬ 
tract properly contributions of higher excited states like the Ai-meson or the 
continuum. 

Note that consideration of the contribution arising from pion exchange 
alone is sufficient as long as we consider only the flavor non-singlet com¬ 
bination f^'^\NS) = Consideration of the flavor singlet 

combination -\- f^‘^\d) naturally forces us to take into ac¬ 

count contributions from r] and g'. In the chiral limit of massless quarks 
which we employ in this paper the g can be considered as massless. The 
axial anomaly, however, will give rise to the g' mass |^. The pedestrian 


9 










solution to this problem would be to disregard the rj' contribution as a short- 
distance one but the real physics is certainly more complicated |^|. As our 
calculation offers no insight into this complicated problem we have decided 
simply to include the rj and the rj' on the same footing as the pion, but it 
should be kept in mind that due to the axial anomaly the prediction for 
is subject to an uncertainty which is probably small, but presently 
unresolvable. 

As far as is concerned an estimate due to Novikov et. al. can 
be found in the literature. We essentially repeated their calculation for 
correlators which determine and with the result = 0.21 GeV^, 
= 0.033 GeV^, = —0.1 GeV^, = 133 MeV. Hence, the contribution 
of a massless chiral boson to the correlation functions at zero momentum can 
be estimated as: 


i J d'^yiOlTOaiy) (g7a75g) (0)|0) = 

i J d'^y{0\TO^{y)i (qgG^^D^^q) (0)|0) = -fl5l5lg^„ 

^ I d%{0\TO^{y) {qa^,g&^D^q) (0)|0)e = -fjldlg^^ (24) 


The factor | in the correlator of with the axial current is due to current 
conservation, see the discussion above. The remaining correlation function 
n = i/(i^|/( 0 |TOo-( 2 /)O'^( 0 )| 0 ) was evaluated with the help of an additional 
sum rule, see |T^, to be H ~ 3 ■ 10“^GeV®. After Fourier transforming 
eq. (1^) and inserting it back into eq. (|15|) we finally obtain the expansion of 
ppB(p 2 ) terms of both local and bilocal power corrections: 


^^(-P^)^log(-pVh^) + {B + Bbl) — <qq>^ (-p^) log(-pVh^) 


TT'^ 


IT 


+C- < fg^GGG > i-p^) log(-p7/i^) 


TT^ 


+DIl^{-p'^) log(-p^/p^) F—ml < qq 


TT^ 


TT 


^G <qq>'<fG^> +iH+ hbb) 

_pZ _pZ 

+l7ias < qq < qq > 


-p- 


p. 


(25) 
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-Bbl 

D 


/ 

J 

K 

d 

-20/ 27 

5/36 

-1/27 

-64/9 

2727 

5/9 

u 

-16/ 27 

1/9 

2/9 

208/9 

2727 

7/36 


Table 2: Numerical coefficients corresponding to BPC contributing to the sum 
rule eq. |2^ . The upper line gives the values for the twist-4 operator involving d 
quarks, the lower line for the corresponding operator with u-quarks 


The numerical coefficients corresponding to the BPC can be read off from 
table H 

Before we proceed to extract the matrix element one comment has to 
be made. An inspection of eq. (pSD reveals an important difference with re¬ 
spect to the previous calculation of |^. The theoretical side of the present 
sum rule given by eq. (^) is manifestly free from effects of mixing of the 


three-point correlation function with two-point correlators [|T^ Q. In 
other words the use of the quark-gluon current o resulted in a much 
milder singularity structure of the three-point correlator than in the case 
with the three-quark current. The mixing, which occurs already at the tree- 
level i.e., before genuine radiative corrections are considered, produces extra 
UV logarithms in Wilson coefficients and makes it necessary to use a more 
complicated model of the spectral representation [^, although the resulting 
corrections turn out to be much smaller than the overall uncertainties. In 
the present case the representation (^) is sufficient to adequately reproduce 
all terms arising from the theoretical calculation. 

To extract the matrix element of interest we employ the standard strategy 
of QCD sum rules. First we multiply the sum rule by m\ — to eliminate 
the single-pole term in (0 and then apply a Borel transformation to both 
sides, arriving at the following expression 


= i(5! - 4! + B{MWj^E2 - 2M^E^) + 

C{M^E2 - M‘^mlEi) + Dml (26) 


where A,B,C and D are the coefficients in front of (—log(— 

(—p^) log(—p^//r^), log(—p^//i^) and l/(—p^) respectively. i?„(so,M^) de¬ 


ll 


























notes as usual 


En = 


n—1 


1 _ ,(-»/«») ^ A . 


fc =0 


(27) 


The dependence on the overlap integrals A/ and Ac can be eliminated from 
eq. (^61) by dividing the three-point sum rule by the two-point sum rule 
derived in ||T6[ 


2(27r)^m^A/AGe-™^/^' = 

5 vr 

+ ^ < g^GG > M^E2 - ^ — (27r)'‘ < qq M‘^Ei + ^(27r)^mo < qq 

2d O TT O 

(28) 

The used set of condensates < qq >= (—0.257 GeV)^, < agj'nGG >= 
0.012 GeV^, ml =< qgaGq > / < qq >= 0.65 GeV^ and < fg^GGG >= 
0.046 GeV® correspond to the standard ITEP values rescaled to the normal¬ 
ization point nl ~ m% ~ 1 GeV^. The strong coupling constant at 1 GeV is 
taken to be as = 0.37 (A = 150) MeV. The continuum threshold is chosen 
as So = (1.5 GeV)^ roughly corresponding to the Roper resonance position. 

The quotient sum rule for the matrix element of the twist-4 operator has 
been plotted in Fig. 2. This figure shows the singlet (S) and nonsinglet (NS) 
part of /(2). /(2)(^) = /(2)(„) + /(2)(iV^) = fi^)(u) - For 

comparison the corresponding sum rules obtained from the analysis in |]^ 
which employed Ioffe currents only are also shown. The square of the overlap 
integral Xj is determined from the additional two-point sum rule 


2(27r)^A?e-”^^/^' = M^E^ 


2 < g^GG > M^Ei + ^{271^ < qq 


(29) 


which is just the standard sum rule considered by Ioffe. Instead of using 
a fixed value for A| we divided the sum rule obtained in [jT^ by the sum 
rule (^). The final values of the matrix elements can be estimated from 
the hgures by taking th m% ~ 1 GeV^ ~ /Xq, where /Xg represents the 
normalization point of the matrix element (®. Numerically, their values are: 


/(2)(^) = _o.09 ± 0.02 f^^\NS) = -0.07 ± 0.02 
f^‘^\proton) = -0.037 ±0.006 
f^‘^\neutron) = —0.013 ± 0.006 


(30) 
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Figure 1: Stability plot of the sum rule eq. (Pq). The full line corresponds to 
the dotted line to For comparison the results of the analysis 

in 1^ are also shown. The dashed line corresponds to fsEK^^) space-dotted 

line to fSU(NSy 
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at /ig 


1 GeV^. These values are to be compared with those obtained from 
the sum rules given in [|T3 : 


/bbk(*^) = -0.068 ± 0.03 fiMNS) = -0.18 ± 0.04 
fiBKiproton) = -0.049 ± 0.01 

/bbk (’neutron) = 0.01±0.01 

All errors given are only due to the dependence on the Borel parameter 
M^. An additional error enters due to the factorization of high dimensional 
condensates. For condensates of dimension 8 the generally accepted error is 
estimated to be of order ~ 20%. There is unfortunately very little experience 
with condensates of dimension 10 since such high dimensions occur seldom 
in calculations. The good agreement of our previous calculation 
results of BBK 


with the 


r3| may be considered as a support for applicability of the 


factorization procedure for dimension-10 condensates as well. In our opinion 
an estimated error of ~ 50% is a very conservative guess. 

Let us now discuss the importance of the various contributions entering 
the expansion (^). As in |T^, Q the sum rule turns out to be dominated 
by the operators of highest dimension i.e., those of dimension 10. To be 
sure that this is a physical effect and not the onset of a breakdown of OPE 
one should reliably estimate the next term in the series, resulting in con¬ 
tributions of dimension 12, which is clearly very difficult. Contrary to the 
evaluation of the twist-3 matrix element where BPC did not contribute, 
local and bilocal corrections enter on equal footing in the present sum rule. 
The numerically important BPC turn out to be those involving the correlator 
with the axial current operator q'yalbQ and the correlator with the operator 
qjbCri^xDaq for which an exact low energy theorem holds. The contributions 
from 6^ and 6^ are numerically smaller and enter with opposite signs so that 
they partially cancel each other. Our hnal results rely on the dimension- 
10 condensates < qq ml, < GG >< qq and the bilocal correction 
< qq i J d‘^y{0\TOfj{y) (q'f'^'y^q) {0)\0). This situation is different from 


^ where the contribution of the dimension-8 condensate < qq ml plays 
the crucial role. It is by no means trivial that a large number of different 
contributions merge together to give a result similar to that of BBK. 

As we have mentioned already the introduction of an explicit gluonic com¬ 
ponent in the nucleon interpolating current resulted in a sum rule which to 
the lowest order is free of extra UV logarithms due to mixing, and therefore 
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the additional uncertainties discussed in |^, |^, Q do not influence the flnal 
estimate. In the previous calculation of the twist-3 matrix element the 
mixing logarithm arose in the dimension-8 contribution which was numeri¬ 
cally negligible in practice. So, as expected beforehand, the consideration of a 
non-perturbative gluonic component has lead to much less singular behaviour 
of the correlators and therefore to more stable numerical predictions. 

In our previous calculation the values of the twist-3 matrix element where 
found to be d^'^\S) = —0.068 ± 0.03 and S) = 0.078 ± 0.03. Using 
these numbers we And from eqn. (^, |^) that both colour electric and colour 
magnetic fields in the rest system of the nucleon contribute at the same order 
of magnitude to the spin 

{pS\gB‘^u'^u\pS) = 

{pS\gB^dU\pS) = 

{pS\ (-U7 u) X gE IpS") = 


{pS\ (dy d) X gE \pS) = 


-(0.07±0.08)m^5° 
(0.188 ±0.08)m^^'" 
(0.09 ± 0.08)m^^'" 
(0.21 ± 0.08)m^5'" . 


(31) 


Obviously such a result shows that simple phenomenological models moti¬ 
vated as analogy to QED are misleading. In any such model one would expect 
the colour-magnetic term to dominate. 

Finally we can analyze the higher-twist contributions to the integral over 
gi{x) and gi~'^{x). If we estimate from experiment, a^‘^\proton) = 
0.022 ± 0.002 from EMC experiments |l| and a^‘^\neutron) = 0.000 ± 0.003 
from E142 we obtain for the Bjorken sum rule 


GeV^ 


(32) 


(33) 


[ dx or"(i, Q“) = la*“>(AfS) + (0.003 ± 0.008) 

Jo 6 

where the twist-2 matrix element is defined as 

2S'’a<")(/) = (re|g,(0)7'’7'<(/(0)|FS) . 

In the case of the difference of proton and neutron structure function the 
flavor index refers to the non-singlet combination u — d. Using isospin 
symmetry this combination can be related to the nucleon /d-decay constant 
QAlgv = a^^'^{NS) leading to the celebrated Bjorken sum rule p. In case of 
the proton spin structure function we get 


I GeV^ 

/ dx gi{x,Q'^) = -a^^\proton) — (0.015 ifc 0.007) 

/o 2 


(34) 
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Here a^^\proton) contains non-singlet and singlet combinations. Using SU(3) 
symmetry it can be expressed by the F and D hyperon decay matrix ele¬ 
ments. To compare with the leading twist matrix elements we take the 
Bjorken sum rule prediction for a^^\NS) and for the Ellis-Jaffe leading twist 
element the measurement of SMC taken at assumably asymptotic . 

dx ^ cx)) = -— = 0.2095 ± 0.0005 

Jo 6 gv 

C dx glix, ~ lOGeV^) = 0.136 ± 0.011 ± 0.011. (35) 

Jo 

Thus we can conclude that the present analysis suggests that the higher-twist 
corrections to both sum rules are small for average in the SMC and EMC 
range 0, §. 
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